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ABSTRACT 

Using simple geometrical arguments, we paint an overview of the variety of 
resonant structures a single planet with moderate eccentricity (e < 0.6) can create 
in a dynamically cold, optically thin dust disk. This overview may serve as a key 
for interpreting images of perturbed debris disks and inferring the dynamical 
properties of the planets responsible for the perturbations. We compare the 
resonant geometries found in the solar system dust cloud with observations of 
dust clouds around Vega, e Eridani and Fomalhaut. 

Subject headings: celestial mechanics — circumstellar matter — interplanetary 
medium — planetary systems — stars: individual (a Lyrae, (3 Pictoris, e Eridani, 
Fomalhaut) 

1. Introduction 



Direct imaging of nearby stars can not yet detect light from extrasolar planets. However, 
imaging can detect circumstellar dust, and when a planet orbits inside a dust cloud, the 
planet can reshape the cloud dynamically, as the Earth perturbs the solar dust cloud. Several 
debris disks around nearby main sequence stars show structures and asymmetries which have 
been ascribed to planetary perturbations (Burrows et al. 1995; Holland et al. 1998; Schneider 
et al. 1999; Koerner et al. 2001; Holland et al. 2003); perhaps these perturbed disks are 
signposts of extrasolar planetary systems. 
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Many of these disk features can be modeled as dust trapped in mean motion resonances 
(MMRs) with a planet. Gold (1975) suggested that as interplanetary dust spirals into the sun 
under the influence of Poynting- Robertson drag (P-R drag), planets could temporarily trap 
the dust in MMRs, creating ring-like density enhancements in the interplanetary cloud. Since 
then, both the InfraRed Astronomical Satellite (IRAS) and the Diffuse InfraRed Background 
Explorer (DIRBE) on the Cosmic Background Explorer (COBE) satellite have provided 
evidence for a ring of dust particles trapped by the Earth (Jackson & Zook 1989; Reach 
1991; Marzari & Vanzani 1994; Dermott et al. 1994; Reach et al. 1995). Models of Kuiper 
Belt dust dynamics (Liou & Zook 1999) suggest that Neptune may also trap dust in first- 
order MMRs. 

Other stars may host planets like the Earth or Neptune. However, most of the known 
extrasolar planets do not resemble the Earth or Neptune; they have masses in the range of 
0.3-15 Jupiter masses, and they often have significant orbital eccentricities (see, e.g., the 
review by Marcy & Butler 2000). Simulations by Kuchner & Holman (2001) show that 
planets as massive as these on eccentric orbits placed in a cloud of inspiraling dust often 
create two concentrations of dust placed asymmetrically with respect to the star. Maps of 
the vicinity of Vega made with the IRAM Plateau de Bure interferometer at 1.3 mm (Wilner 
et al. 2002) and with the JCMT at 850 /xm (Holland et al. 1998) reveal two concentrations 
of circumstellar emission whose asymmetries can be naturally explained by such a model, 
possibly indicating the presence of a few- Jupiter mass planet in an eccentric orbit around 
Vega (Wilner et al. 2002). Other papers have numerically explored the interactions of par- 
ticular planetary system configurations with a dust disk, with a view towards developing a 
general key for interpreting disk structures (Roques et al. 1994; Lecavelier Des Etangs et al. 
1996; Liou & Zook 1999; Quillen & Thorndike 2002). 

We assemble a primitive version of such a key by mapping the geometries of the MMRs 
which are likely to trap the most dust near a planet embedded in a debris disk. We illustrate 
the patterns formed by the libration centers of the trapped particles in an inertial frame — the 
frame of a distant observer. These basic patterns allow us to characterize four structures 
that probably span the range of high-contrast resonant structures a planet on an orbit with 
eccentricity up to ~ 0.6 and low inclination can create in a dust disk. 

Figure 1 shows these structures here for reference; we discuss them throughout the paper, 
particularly in Section 4. Cases I and II in this figure represent the structures formed by 
planets on low-eccentricity orbits. Cases III and IV represent structures created by planets 
on moderately eccentric orbits. Cases I and III represent structures created by planets with 
substantially less than 0.1% of the mass of the star; cases II and IV represent structures 
created by more massive planets. 
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Fig. 1. — Four basic resonant structures: I) low mass planet on a low eccentricity orbit, 
II) high mass planet on a low eccentricity orbit, III) low mass planet, moderate eccentricity 
orbit, and IV) high mass planet, moderate eccentricity orbit. 
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2. 



PLANETS ON LOW-ECCENTRICITY ORBITS 



2.1. 



Low Mass Planets 



First, we review some of the physics of resonant dust rings created by relatively low 
mass planets, like the Earth and Neptune. In the rest of this paper, we generalize the 
discussion to planets with mass up to ~ 15 Jupiter masses and eo up to ~ 0.6. The resonant 
geometries we will describe apply to the general restricted three-body problem. However, 
we have in mind a a source of dust, like the asteroid belt or the Kuiper belt, which releases 
dust on low-eccentricity orbits at semimajor axes, a > do, where a is the semimajor axis of 
the particle's orbit and do is the semimajor axis of the planet's orbit. We use the convention 
that quantities with a subscript refer to the planet. 

We do not consider interactions between the dust grains, or between dust and any gas 
in the disk. For example, our approach does not apply to the gaseous, optically-thick disks 
found around young stellar objects. However, the resonances we describe must underlie 
the basic resonant features of the solar system dust complex and debris disks around main 
sequence stars with less than a few lunar masses of dust. 

The radiation forces on a particle are parametrized by /3, the strength of the stellar 
radiation pressure force on a particle divided by the strength of the stellar gravitational 
force on a particle. For a spherical particle with radius s > 1 fim and density 2 g cm~ 3 
orbiting a star with mass M* and luminosity L*, 



A typical particle in the solar system may have s ~ 1-100 (Grun et al. 1985; Fixsen & 
Dwek 2002), or j3 ^0.285-0.00285. Dust grains released at circumstellar distance r that are 
too large to be ejected by radiation pressure (j3 < 0.5) spiral into the star via P-R drag 
(Robertson 1937; Burns et al. 1979), on a time scale 



Other drag forces, like solar wind, may also contribute to the decay of the particle's orbit 
(Banaszkiewicz et al. 1994). 

Dust spiraling inwards towards a planet encounters a series of exterior MMRs, each of 
which is associated with terms in the disturbing function of the form 




(1) 




(2) 
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F(a, e, e ) cos0, 



(3) 



a 



- 5 - 



where G is the gravitational constant, m is the mass of the planet, and a = a^/a (see, e.g., 
Brouwer & Clemence 1961; Murray & Dermott 1999). The resonant argument, 0, is a linear 
combination of the orbital elements of the particle and the planet, which can be interpreted 
as an angle (Greenberg 1978). We do not discuss interior MMRs because they can not sustain 
long-term trapping (see e.g. Murray & Dermott 1999, p. 381). The potential, < R > re s, 
causes the argument to accelerate. For a particle trapped in the resonance, librates 
about pa if F(a, e, e ) < 0, and librates about pa -k if F(a, e, e ) > 0. A MMR of the 
form j:k is nominally located at a semimajor axis given by l/a pa (j//c) 2//3 (l — /3) 1 / 3 . 

As a particle approaches the planet from afar, it encounters stronger and stronger reso- 
nances, and has a better and better chance of becoming trapped. For particles approaching 
the Earth or Neptune, the first resonances which are strong enough to trap substantial 
amounts of dust are the first-order MMRs: resonances of the form 2:1, 3:2, 4:3, etc. To 
first order in eccentricity and inclination, these resonances consist of pairs of terms with 
arguments 

01 = J'A - (j - l)Ao - C7 (4a) 

02 = J'A - (j - l)Ao - zu , (4b) 

where A and Ao are the mean longitudes of the particle and the planet, and w and wq 
are the longitudes of pericenter of the particle and the planet, respectively. For the Earth 
(e = 0.017) and Neptune (e = 0.0086), the 0i resonance dominates at all values of zu 
if the dust particle has even a small orbital eccentricity, and 0i librates about 0i ~ n. 
Passage through this resonance slowly raises the particle's eccentricity, which asymptotically 
approaches a limiting value, e max , when the dynamics are followed using an expansion to 
second order in e (Weidenschilling & Jackson 1993; Sicardy et al. 1993; Beauge & Ferraz- 
Mello 1994; Liou & Zook 1997): 

e max = y/2/(5j). (5) 

In the process, the particle's orbit becomes planet crossing, and in a matter of a few P- 
R times, the particle generally leaves the resonance via a close encounter with the planet 
(Marzari & Vanzani 1994) and an abrupt transition to a new orbit with a different eccentricity 
and semimajor axis. Near the Earth, trapped dust populates several first-order MMRs. 

The condition that 0i librates about 0i ~ it provides a relationship between A and w. 
Since rotations of the whole system by 2ir must not affect the dynamics, we can write that 
the libration centers are located where 

A pa ((j - 1)(A + 2-kK) +w + for K e Z. (6) 

Figure 2 shows how this condition leads to the formation of a density wave, using the 3:2 
MMR as an example. Figure 2a shows a variety of elliptical dust orbits, all with the same 
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1/a = 1.31 and e = 0.8e max , but with different ca's. According to Equation 6, each one 
of these orbits has j '• = 3 different longitudes that may be libration centers for a particle 
trapped in the <f>i term. Figure 2b shows these three locations, marked with X's, for each of 
two elliptical orbits with slightly different m's. 

Figure 2c shows the locus of the libration centers for the whole range of elliptical orbits 
shown in Figure 2a. In Figure 2d, the planet has moved to a new longitude, and following 
Equation 6, so have the libration centers. The density wave formed by the trapped particles 
resembles the pattern formed by the libration centers, blurred somewhat by libration, mostly 
in the azimuthal direction. Ozernoy et al. (2000) shows some examples of how libration blurs 
the patterns formed by particles trapped in individual resonances. 

Though any individual particle in Figure 2 moves steadily counterclockwise, slower than 
the planet, the density wave appears to rotate as a fixed pattern together with the planet. At 
conjunction (A = Ao) the resonance condition implies that w ~ Ao+tt, so the libration centers 
just outside the planet are always near apocenter. This condition creates the signature gap 
at the location of the planet seen in simulations of the Earth's ring (Dermott et al. 1994) 
and of Kuiper Belt dust interacting with Neptune (Liou & Zook 1999) and illustrated in 
Figure 1 (Case I). 

2.2. Higher Mass Planets 

We can understand the range of possible resonant dust cloud structures by understand- 
ing the density-wave patterns created by the series of MMRs that are likely to trap dust 
near a planet. Table 1 lists the arguments in a series of MMRs that figure most prominently 
in the sculpting of dust clouds (ignoring the planet's inclination, which first appears with 
order inclination squared) in order by a. It also lists the corresponding leading terms in 
F(a, e, eo), evaluated for \x = 1/1047 (i.e. Jupiter), and f3 « 1. The first MMRs in the list 
are the first-order resonances described above. 

For a given planet/dust cloud combination, a few MMRs generally dominate the ob- 
served resonant structure. In simulations of /3 — 0.037 dust particles approaching the Earth 
(Dermott et al. 1994), the first-order j =4-10 MMRs dominate the appearance of the trapped 
dust cloud. In simulations of Neptune's ring by Liou & Zook (1999), the 2 : 1 and 3 : 2 
resonances (j = 2 and j = 3) dominate the appearance of the dust cloud. 

The factor of Gmo/a in Equation 3 makes MMRs at large a weaker and less able to trap 
substantial quantities of dust. But for more massive planets, the trapping probability for all 
MMRs is higher (Lazzaro et al. 1994), so particles approaching from afar become trapped 
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Fig. 2. — How a planet on a low eccentricity orbit creates density waves, a) Several particle 
orbits with different re's, b) Libration centers of the 3X — 2X — w term on two of these orbits, 
indicated by X's. c) The locus of all the libration centers, d) The density wave appears to 
orbit at the same angular frequency as the planet (located at the filled circle). 
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sooner, in MMRs with longer orbital periods and larger semimajor axes. Likewise, close 
encounters with the planet more easily scatter dust grains out of exterior MMRs at small a, 
and this effect worsens with the mass of the planet. 

The MMRs near the planet also become close together and begin to compete with one 
another. If \x = m /M±, the resonance overlap criterion (Wisdom 1980; Duncan et al. 1989) 
predicts that first-order resonances with j > 0.45 (/x(l — f3))^ 2 ^ 7 + 1 are completely chaotic. 
This condition appears to place an absolute limit on how far down the chain of resonances a 
dust particle can be trapped. For large particles near the Earth (// pa 3 x 10~ 6 ), this criterion 
predicts that the first completely overlapped MMR is 17:18, for Neptune (/i ~ 5 x 10~ 5 ), it 
is 8:9, and for Jupiter (// pa 10~ 3 ) it is 4:5. 

As we have mentioned, many of the observed extrasolar planets have substantially more 
mass than planets in the solar system. Such massive planets quickly scatter dust from 
their first-order MMRs. So after the first-order resonances, Table 1 lists the n:l resonances, 
the lowest order terms available at large semimajor axes — beyond the 2:1 MMR (nominally 
located at a/a = 1.59). In numerical integrations, Kuchner & Holman (2001) and Wilner 
et al. (2002) found that dust spiraling inwards towards a massive planet became trapped in 
this series of n:l resonances. Table 1 shows that the terms in these resonances which don't 
depend on the eccentricity of the planet have the form jX — A — (j — l)w. These terms must 
dominate when eo is small. 

The jX — A — (j — l)w terms in Table 1 are mostly positive, so most of these arguments 
librate around n. However these terms all include contributions from indirect terms, which 
arise from the reflex motion of the star. The indirect contributions tend to reduce the 
strength of the resonance, and in the case of the 3:1, this effect makes the coefficient negative, 
so the 3A — A — 2w argument librates around 0. 

The libration centers for the terms which librate around it are located at 

Aw ((j -l)w + X + ir(l + 2K))/j, for K e Z. (7) 

Again there are j libration centers spaced evenly around each elliptical dust orbit, and again, 
the locus of libration centers appears to co-rotate with the planet. But this time, when we 
set A = Ao, we find that at conjunction, the particle can be in one of a few different places; 

zu pa A — 7r(l + 2K)/(j — 1) at conjunction. (8) 

One can also show following Weidenschilling & Jackson (1993), for example, that passage 
through a pure jX — Ao — (j — l)w term raises the eccentricity of a dust particle in the same 
way that passage through a j X — (j — 1)A — vo term does, and that the limiting eccentricity 
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is again given by Equation 5. We can generate a rough picture of the density waves created 
by the dust trapped in these terms in the same manner as Figure 2, by drawing a variety 
of elliptical orbits with eccentricity e ~ e max at an appropriate semimajor axis, and using 
Equation 7 (and Kepler's equation) to locate the libration centers on these orbits. 

Figure 3 shows the patterns created by the most important resonant terms listed in 
Table 1. The first column of Figure 3 shows patterns for the terms which appear when the 
planet's eccentricity is low: the jX — (j — 1)A — w terms for first-order resonances, and 
the jX — Xq — (j — l)w terms for n:l resonances. The patterns are the locii of the libration 
centers, generated in the manner of Figure 2. 



The patterns in the first column of Figure 3 appear in the textbook by Murray and 
Dermott (1999, p. 325), derived in a slightly different context — by tracing the path of a 
particle in a frame corotating with the planet. However, the literature offers little discussion 
of resonant structures created by planets on eccentric orbits, and this case is crucial for 
understanding extrasolar planetary systems. Naturally, for a planet on an eccentric orbit, 
there is no simple corotating frame. But we can deduce the resonance patterns associated 
with planets on eccentric orbits by building on the arguments used above. 



Besides resonant perturbations, planets introduce secular perturbations to the orbits of 
nearby particles. When the planet's orbit is eccentric, these perturbations generally produce 
a correlation between e and w. We must consider these effects when we examine MMRs 
with planets on eccentric orbits. 

We can demonstrate the importance of secular perturbations by comparing them with 
radiation forces. Near a planet, the P-R time, and thereby the trapping time, is longer than 
the secular time scale for dust grains with (5 < fi , where 



3. 



PLANETS ON ECCENTRIC ORBITS 



3.1. 



Secular Effects 





and the Laplace coefficient (size ~ unity) is 
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Fig. 3. — Patterns formed by dust in MMRs listed in Table 1. 
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I.e., (3 ~ 100/i. For example, /3 = 0.022 for Neptune, (3 = 0.44 for a Jupiter-mass planet 
orbiting Vega at a = 40 AU. So we can expect most observable particles in the cases we are 
interested in (10~ 5 < < 10~ 2 ) to suffer significant secular evolution while they are trapped 
in a MMR. Secular perturbations from planets with eccentric orbits affect all dust particles 
in their vicinity, even those which are not in strong MMRs. 

In the Laplace-Lagrange description of secular perturbations, valid for e, e << 1, a 
particle's osculating e and w are expressed as a combination of a constant forced elements, 
e forced and Wf orC ed, and time- varying free elements, e/ ree and Wf ree . Secular evolution is 
easily visualized in the (h, k) coordinate system, where h = e cos w and k = e sin w. In this 
system, we write h = h forced + h free , and k = k forced + k free . 

When there is only a single perturbing planet on a fixed orbit, the particle's forced 
elements are constant, and the osculating elements, h and k, trace out a circle centered 
on {hf orced , kforced) with a radius of c free- The forced longitude of pericenter is w forced — 
Skictoxi{k forced/ hforced) = ^o, and the forced eccentricity is 



As a approaches a , ef orC ed approaches e . Table 1 lists the approximate value of e/ orce d/e 
at the nominal a for each MMR. 

In a cloud of many particles, many orbits with a range of hf ree and kf ree are occupied. 
However, all particles with a given semimajor axis will have the same h forced and k forced- 
For example, if all particles are released on circular orbits outside the planet's orbit and 
outside of any MMRs, then the Laplace-Lagrange solution prescribes that at any given time 
all the particles with a given semimajor axis will occupy a circle with radius e/ ree centered 
on the point (hf or ced, kf or ced) in the (h, k) plane. Dermott et al. (1985) and Wyatt et al. 
(1999) showed that a set of orbits occupying such a circle in the (h, k) plane form a cloud 
that is roughly circular, but the center of the cloud is offset from the star a distance eocto 
along the planet's apsidal line in the direction of the center of the planet's orbit. Hence, the 
background dust cloud in the vicinity of a large planet should often appear circular in the 
absence of MMRs, though if the planet's orbit is eccentric, the center of the circle will be 
offset from the star. 

When the perturber has higher eccentricity, or when a particle is in a MMR, the particle's 
secular trajectory changes somewhat; the orbit no longer traces an exact circle in the (h, k) 
plane, even when the librations are averaged away (Wisdom 1983). However, the character 
of the averaged secular motion often remains the same; h and k follow a simple closed loop 
around {hforced-, kforced)- We will retain the spirit of the Laplace-Lagrange approximation for 




3/2 
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the secular motion of the particles for the remainder of this paper, and we wil refer loosely 
to free and forced elements, even when we are discussing resonant orbits. We find that this 
approximation suffices for our broad exploration of planetary signatures in debris disks. 



3.2. Mean Motion Resonances 



When a particle is trapped in a MMR with a planet on a circular orbit, the resonant 
argument which only depends on zu, not on zu dominates the particle's motion. The other 
terms in the resonance have coefficients F(e, eo) that are zero when eo is zero. However, when 
a particle is trapped in a MMR with a planet on an eccentric orbit, resonant arguments which 
depend on zu may come into play. 

To leading order, the averaged disturbing function for a MMR with a planet on an 
eccentric orbit is the sum of two or more terms: 



< R>, 



GlTLr 



y^Fg(o;,e,eo) cos0 f , 



where, as Table 1 shows, (pi = p\ — q\o + (p — q)w and 
Wisdom (1983), we can re-express this sum as 



(12) 

i + (£ — l)(zu — voq). Following 



GrriQ 



F'(a, e, e , zu, zu ) cose//, 
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where 



and 



4>' = 0i + arctan 



Eg = 2l^sin((e-l)(P7-s7 ) + *g) 
1^1 + E? =2 1^1 cos ((£ -l)(w- tu ) + St) 



F'(a, e, e , zu, zu ) 



^F f (a,e,e ) sin0 5 



^F ? («,e,e ) cos 5 



(14) 



(15) 



The quantity 5^ = vr for F$ < 0, and <5g = otherwise. We define F' to be always > 0, so on 
resonance, the new argument, 0', librates around 0' ~ ir. Equation 14 shows explicitly that 
the differences among the terms only appear on a secular time scale, as zu — zuo precesses. 

At any moment, the particle's orbit may be viewed as undergoing libration about 0' = n. 
This change in the effective resonant argument can change the constraints on the orbital 
elements of resonant objects, which can result in dramatically different-looking clouds of 
trapped particles. Often one term dominates, in the sense that 0' 0^ or 0' ps 0^ + tc. This 
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circumstance depends on e, eo, to, and w , and the details of trapping, but we can appeal to 
numerical simulations and use what we know about the secular evolution of dust particles 
to help decide which terms are likely to be most important. 

3.3. High Mass Planets 

It is easier to consider high mass planets first, because the n:l resonances that dominate 
the appearance of dust clouds containing a massive planet are farther from the planet than 
the first-order resonances which dominate in the case of a low mass planet. Particles released 
on circular orbits approach a planet with e/ ree ~ ef orce d, so they reach the n:l MMRs with 
e ranging from to ef orce d + e/ ree . As Table 1 shows, ef orce( i ~ e /2 near the lowest order 
n : 1 MMRs. When w pa zuq + 7r, all the terms become degenerate, in the sense that all 
of their libration centers occur at the same longitudes. At the other end of the range of 
secular motion, w ~ zu , and e ~ e , so according to Equation 14, the term with the largest 
coefficient in F(a, e, e ) will dominate the motion of the particles. 

The arguments which match the terms with the largest coefficients in F(a, e, eo) (see 
Table 1) share a second common form: jX — Xq — wq — (j — 2)w. Numerical integrations 
(Wilner et al. 2002) confirm that these terms to dominate the resonant dynamics over a wide 
range of planet eccentricities. These terms always have F(a, e, e ) < 0, so the corresponding 
arguments librate about 0. The j libration centers on each ellipse are located at 

A pa (Xq + wq + (j — 2)w + 2nK)/j, for K e Z. (16) 

We reserve a discussion of the complicated secular dynamics of these resonances for a future 
paper. 

Figure 4 illustrates the pattern formed by particles in the 3X — X — w — w term 
near a planet with e = 0.6. Once a particle is trapped in resonance, its secular evolution 
no longer obeys the Laplace-Lagrange solution, but the libration centers of the trapped 
particles will still occupy a closed curve in the (h, k) plane. Figure 4a shows such a variety 
of orbits centered on h = 0.58eo, k — 0. Orbits which are close to apse-alignment with the 
planet have higher eccentricities. Like the Laplace-Lagrange orbit distributions discussed by 
Dermott et al. (1985) and Wyatt et al. (1999), these orbits taken together form a pattern 
which is roughly circular, but the center of the circle is offset from the star along the planet's 
apsidal line in the direction of the center of the planet's orbit. 

Figure 4b shows the 3 libration centers on two of these orbits, calculated from Equa- 
tion 16. By definition, A = M + w, where M is the particle's mean anomaly. For particles 
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Fig. 4. — Libration centers of the 3A — Ao — vjq — w term, a) Several particle orbits with 
different e and w. b) The libration centers on two of these orbits when the planet is at 
pericenter. c) All the libration centers, d) The clumps formed by particles trapped in this 
term appear to rotate at half the angular frequency of the planet. 
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at pericenter, M = and A = zu, so for particles trapped in any jX — \ — zo — (j — 2)w 
term, 

w « M /2 + zu + ttK iorKeZ, (17) 

where M is the planet's mean anomaly. In other words, the libration centers for this family 
of terms reach pericenter at two different longitudes, and these special longitudes precess at 
an angular frequency equal to half the Keplerian angular frequency of the planet. When the 
planet reaches pericenter (M = 0), so do the particles on orbits that are apse-aligned or 
anti- apse-aligned with the planet's orbit. Figures 4c and d show that the locus of libration 
centers makes a characteristic two-lobed pattern, which appears to rotate at half the mean 
angular frequency of the planet. The center of this apparent rotation lies on the planet's 
apsidal line between the star and the center of the planet's orbit. 



3.4. Low Mass Planets 

The resonances populated near low-mass planets cluster at semimajor axes near the 
planet's semimajor axis, where ef orce d ~ e . In this vicinity, the resonant orbits can easily 
become planet-crossing unless they are roughly apse-aligned with the planet's orbit, i.e. 
^free ^ & forced] and ZU ~ w . Figure 5a shows a collection of roughly apse-aligned orbits, 
all with ef orced = 0.5, and e/ ree = 0.12. In the absence of MMR's, the combination of these 
elliptical orbits would appear as an elliptical ring. 

For the case of the 3:2 MMR, a single one of these orbits has three libration centers, 
shown in Figure 5b. For a family of orbits that must remain roughly apse-aligned with 
the planet, the secular spread in the libration centers is necessarily small, and the resonant 
arguments are all roughly equal. So, as Figure 5c shows, we should expect the density wave 
to resemble the pattern produced by the libration centers in any one orbit; the libration 
centers associated with any MMR of the form j:k form a set of j clumps. As Figure 5d 
shows, the locus of all the libration centers orbits at k/j times the angular frequency of the 
planet; when the planet's mean anomaly changes by tt/2, the mean anomalies of the clumps 
of particles in the 3:2 MMR change by n/3. 

An exterior first-order MMR has two relevant terms (Equations 4a and 4b), and for the 
second term, F 2 (a,e,eo) < 0. For this case, Equation 15 becomes 



0' = 0i + arctan 



\F 2 (a, e, e )| sin(cu - zu + n) 



\F ± (a, e,e )| + [F 2 (a,e,e )\ cos(w - w + n) 



A key difference between the two terms is that one argument librates around 0, the other 
around ir. 
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Fig. 5. — How dust confined to the 3:2 MMR near a planet on a moderately eccentric orbit 
might appear, a) Several orbits with different e and zu, distributed in a Laplace-Lagrange 
circle with small e/ ree . b) The libration centers on two of these orbits when the planet is at 
pericenter. c) All the libration centers, d) The three dust clumps really do orbit at 2/3 of 
the planet's orbital frequency. 
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Deciding which resonances dominate in the vicinity of a low mass planet on an eccentric 
orbit can be tricky. When |F 1 (a, e, e )| >> \F 2 (a,e,e )\, fa ~ fa. When |Fi(o!, e, e )| << 
|F 2 (a;, e, e )|, fa ~ fa + tt- However, when Fi(a,e,eo) ~ F 2 (a,e, e ), the two terms tend to 
cancel each other, that is, \F'\ m \Fi — F 2 \ when w m zu . This cancellation diminishes the 
strength of the resonance. 

So we should expect first-order resonances to be relatively weak for a dust particle 
whose orbit has been secularly apse-aligned with the planet's orbit. At e = ef orce d, the fa 
resonance generally has the larger coefficient. The 2:1 MMR is an exception, because i*\ 
for the 2:1 is diminished by an indirect term, so the fa resonance easily dominates. We 
leave it to numerical studies (Quillen & Thorndike 2002) to decide which MMRs and terms 
dominate in a given system containing a low mass-high eccentricity planet. However, no 
matter which MMRs dominate, we expect the same generic behavior; the trapped dust will 
form an eccentric ring of dust clumps. 



4. PLANETARY SIGNATURES 

4.1. Four Basic Cloud Structures 

The first column of Figure 3 illustrates the density waves formed by the resonant terms 
we expect will dominate the appearance of dust clouds near planets with low orbital ec- 
centricity. The second and third columns illustrate the density waves formed by the res- 
onant terms we expect will dominate near planets with moderate orbital eccentricity: the 
jX — (j — l)Ao — vd§ terms for first-order resonances, and the jX — Xq — voq — (j — T)w terms 
for n:l resonances. In the second column, the planet is at pericenter; in the third column, 
the planet is at M = tt/2. 

Our map of these basic patterns prepares us to consider more generally the appearance 
of a dust disk in the vicinity of a planet, where trapped dust occupies several MMRs. Now 
we can describe the origin of the structures shown in Figure 1; they are superpositions of 
the patterns shown in Figure 3. 

Cases I and II are superpositions of the patterns depicted in the left-hand column of 
Figure 3. Cases III and IV are superpositions of the patterns depicted in the right-hand 
column of Figure 3 (and also in the middle column, at a different planetary orbital phase). 
In general, the rings created by the low-eccentricity planets appear to co-rotate with the 
planet, but the resonant structures created by moderate-eccentricity planets do not, because 
the density wave patterns associated with the resonances they excite vary with the planet's 
orbital phase. 
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Case I: A low mass planet with low orbital eccentricity, like the Earth or Neptune, traps 
dust in first-order jX — (j — 1)A — w resonances. Case I in shows a superposition of patterns 
produced by terms of this form, the patterns in the upper left-hand column of Figure 3. 
Populating these resonances creates a ring with a gap at the location of the planet. 

Case II: A higher mass planet on a low-eccentricity orbit traps dust in more distant n:l 
resonances, in terms of the form jX — A — (j — l)zu. Case II shows superpositions of patterns 
produced by such terms, the patterns in the lower left-hand column of Figure 3. These 
resonances create a larger ring with a smooth central hole. 

Case III: A low mass planet on a moderately eccentric orbit traps dust in MMRs with small 
secular motion near apse- alignment with the planet, creating a blobby eccentric ring. The 
blobs in this ring appear to be continually created and destroyed, as dust clumps occupying 
different MMRs pass through one another. The more highly concentrated clumps are located 
near the apocenter of the planet's orbit, where all the particles in the jX — (j — l)Ao — t^o 
terms (except for the 2:1) are constrained to have their conjunctions. The planet may or 
may not be located near a gap in the ring. 

Case IV: A high mass planet on a moderately eccentric orbit creates a ring offset from the 
star containing a pair of clumps, where the two-lobed patterns of all the n:l resonances of 
the form jX — Ao — (j — 2)w — zu conincide. The clumps appear to orbit a point along the 
planet's apsidal line, between the star and the planet, at half the mean orbital frequency of 
the planet. 

Naturally, intermediate cases will result in more than one variety of resonance being 
populated by dust particles; these clouds can possess features of more than one of the extreme 
cases shown in Figure 1. For example, an intermediate mass planet on a low eccentricity 
orbit will have a slightly larger ring with a less prominent gap at the location of the planet 
than Case I. A planet with intermediate eccentricity may have both a ring component and 
blobs, though librational motion may smear these blobs so that they blend into the ring. 



4.2. Strong Drag Forces; Wakes 



When the drag force is large and the planet's mass is small, the libration centers can 
shift substantially from ir. For a MMR containing a collection of resonant arguments of the 
form (f) — jX — kXo — (j — l)zu — — 1)(zuq — zu), the resonant perturbations are 

da 



dt 



—2janfiF'(a, e, e , w, w Q ) sin< 



(19) 



Res 
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de 
~dl 



—nfiF'(a, e, e , w, zu ) sin 0', 



Res 



where the particle's angular frequency, n, on resonance is nominally 

n=(GM*(l-P)/a 3 ) 1/2 . 
PR drag causes relatively slow changes in a and e (Wyatt & Whipple 1950): 



da 
~dl 
de 
dt 



PR 



PR 



-GM+(3 (2 + 3e 2 ) 
a~c (1 - e 2 ) 3 / 2 
-5GM*(3 (2 + 3e 2 ) 



2a 2 c 



e 2)l/2 



(20) 

(21) 

(22) 
(23) 



where c is the speed of light. To locate the libration centers, we set [da/dt] PR + [da/ dt] PR = 0, 
and find 

. , (GMjay/i(3(l-{3) (2 + 3e 2 ) 

8111 <Po " 2j/,cF'(a, e, e , U7 ) (1 - e 2 ) 3 / 2 ' 1 J 

On resonance, 0' librates around O , which is greater than ir. 

For the case of a planet on a circular orbit, the observable effect on a dust cloud is a 
shift in the locations of the pericenters of the orbits of the trapped particles. In the absence 
of this effect, the libration centers first reach pericenter an angle n/k behind the planet. 
With this effect, the first pericenter is located closer to the planet, at an angle (2ix — <po)/k. 
This asymmetry concentrates the trapped dust from the Earth's ring into a blob trailing the 
Earth, sometimes called the Earth's "wake" (Dermott et al. 1994; Reach et al. 1995). 

Figure 6 shows an example of how cases I and III, as illustrated in Figure 1, could 
appear when the perturber's gravity is weak compared to the drag force and the shifts in 
the libration centers are substantial. Case I shows a wake trailing the planet. For the terms 
illustrated in columns two and three of Figure 3, the shift in the libration center appears 
as a displacement of the density enhancements. The n:l resonance clumps will appear at 
w = Kir + Mo/2 + wq + (0o — 7r )/2, i.e. their longitudes shift by (0o — 7r )/2- The first-order 
resonance clumps near low-mass planets on moderately eccentric orbits (Case III) shift by 
(0o — Tr)/k in the prograde direction. 

This effect depends on the velocity of the dust (~ (GM^/a) 1 / 2 ), on /3, and on the mass of 
the planet; it is more pronounced for small dust grains and low-mass planets close to massive 
and luminous stars. It is negligible for patterns generated by more massive planets (cases 
II and IV) , even if the planets are in 3-day orbits. Moreover, this effect does not appear in 
simulations of Neptune's ring (Liou & Zook 1999) since even Neptune has high mass (and 
low orbital velocity) compared to the Earth. The e Eridani ring seen by the JCMT (Holland 
et al. 1998) probably does not show this effect either for the same reasons. 



-20- 




Fig. 6. — Cases I and III when the drag force is particularly strong and the planet's mass 
is small. Case I develops a "wake" trailing the planet. The blobs in Case III shift in the 
prograde direction. 
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5. DISCUSSION 

5.1. Observed Examples 

Earth is the only planet in the solar system whose resonant ring has been detected 
(Kuchner et al. 2000). This wake in this ring produces an asymmetry in the thermal infrared 
background between the direction leading the Earth and the direction trailing the Earth, 
which IRAS and COBE DIRBE detected (Reach 1991; Dermott et al. 1994; Reach et al. 
1995). Indeed, the Earth's ring is the only observed resonant dust cloud structure — in 
the solar system or elsewhere — identified with an independently observed planet! However, 
models of Kuiper belt dust interacting with the outer planets suggest that Neptune also 
creates a type I ring with a characteristic hole at the location of the planet. 

Images of nearby stars with debris disks supply further examples of rings which point to 
undetected planets. Dust clouds around Vega, j3 Pictoris, e Eridani and Fomalhaut have all 
been imaged at submillimeter wavelengths (Holland et al. 1998; Greaves et al. 1998; Holland 
et al. 2003). The (3 Pictoris disk is edge-on, making its resonant rings hard to distinguish — if 
it has any. However, the other systems appear to contain resonant rings which may be easier 
to classify. 

Vega: The Vega debris disk may provide an example of case IV. Submillimeter images of 
Vega (Wilner et al. 2002) show two concentrations of emission that resemble those in Figure 1; 
they are not co-linear with the star, and one is closer to the star than the other. Wilner et 
al. (2002) used this model to estimate that the perturbing planet has mass fi ~ 10~ 3 and 
eccentricity e ~ 0.6, and used the locations of the two knots of emission to infer the planet's 
current loation. 

e Eridani: Case III may be a good model for the e Eridani dust ring, which has roughly 
four major concentrations of emission irregularly placed around the ring, and an apparent 
gap in the ring. As Figure 1 shows, the low-mass, moderate eccentricity planet model 
naturally explains the presence of a few dust clumps of various concentrations. Quillen & 
Thorndike (2002) have developed a numerical model of type III for the e Eridani dust ring 
using the numerical technique developed by Kuchner & Holman (2001). They find that dust 
released just outside a planet with eccentricity 0.3 and mass /i = 10~ 4 becomes trapped in 
the 5:3 and the 3:2 MMRs with the planet, and forms a blobby eccentric ring like the one 
in Figure 1. The second order 5:3 MMR may be relatively strong because the first order 
resonances suffer from the cancellation described in Section 3.4. Though dust trapped in 
these resonances probably forms relatively narrow rings, the large beam of the JCMT could 
easily blur such a narrow blobby ring into something resembling the Holland et al. (1998) 
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image, as Quillen & Thorndike (2002) illustrate. 

Fomalhaut: The A3 V star Fomalhaut has a circumstellar dust cloud, recently identified as 
ring with an azimuthal asymmetry representing ~ 5% of the total flux from the disk (Holland 
et al. 2003). This ring is difficult to classify, since it is roughly 20° from edge-on. Wyatt & 
Dent (2002) suggest that the dust may trapped in a 2:1 MMR with a planet, and Holland 
et al. (2003) suggest that the dust might be trapped in a 1:1 MMR with a planet. However, 
the 1:1 MMRs are much harder to populate with dust than the MMRs we consider here, and 
contrary to Ozernoy et al. (2000), we consider it improbable that a single MMR, even the 
2:1, could dominate the appearance of a real dust cloud, particularly one as collisional as 
Fomalhaut's. More likely, the structure in the Fomalhaut ring arises from collisions Wyatt 
& Dent (2002) or if resonant trapping is responsible for the azimuthal structure, there is 
another clump to be found, perhaps obscured by the limb-brightening. 

We know no observed examples of case II rings. However, the close-in extrasolar planets 
with periods less than ~ 30 days have nearly circular orbits, probably due to tidal effects 
(Marcy & Butler 2000). Any dust disks associated with these planets could form examples 
of case II. 



5.2. Other Considerations 

We have only addressed the geometry of resonant dust orbits. A variety of other phe- 
nomena may affect the appearance of an actual debris disk. For example, we have not 
discussed how the resonances are populated. Different resonances may dominate when a 
disk is fed by a belt of asteroids or comets whose orbits are restricted to a small region of 
phase space. 

Furthermore, many known dust disks are collisional. For systems like the rings of Saturn, 
where each particle has a collision roughly once per orbit, the collisions dominate the resonant 
effects we catalog here, and a fluid description of the particles becomes more appropriate. 
Collisions which occur on intermediate time scales, however, may leave particles trammelled 
in the net of the underlying strong resonances, altering only the spectrum of populated 
MMRs. For example, mutual collisions among dust grains may destroy dust particles before 
they can access all the strong resonances; this effect would preserve the basic four structures 
shown in Figure 1. 

We have also neglected the terms in MMRs that depend on the inclinations, and re- 
stricted our analysis to planets with moderate eccentricity. At high planet inclinations 
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(i 2 ~ 1) and eccentricities (e 2 ~ 1), many new terms in the disturbing function become 
relevant. These effects can alter the basic resonant geometries. 

Many planetary systems have more than one planet; we have also only considered the 
effects of one planet. Since we are only interested in high-contrast structures, we may 
justify this approach by saying we only care about the first massive planet to encounter the 
inspiraling dust — massive enough to create a high-contrast structure by ejecting most of the 
dust grains as they pass. For example, in simulations by Liou & Zook (1999), Neptune ejects 
most of the inspiraling Kuiper belt particles before they can encounter any other planets; 
we are interested in analogous planets. However, the secular evolution of a multiple-planet 
system may affect how a planet interacts with a dust cloud, even if little dust reaches most 
of the planets, and they do not carve their own signatures. 

Finally, we have tacitly assumed that the dust cloud is observed face on. Disks which 
are tilted from face-on may show a variety of asymmetries due to effects other than resonant 
trapping, such as the IRAS/DIRBE dust bands (Hauser et al. 1984), secular warps (Augereau 
et al. 2001), and limb brightening. 

6. CONCLUSIONS 

Four basic structures probably represent the range of high-contrast resonant structures 
a planet with eccentricity < 0.6 can create in disk of dust released on low-eccentricity orbits: 
a ring with a gap at the location of the planet, a smooth ring, a blobby eccentric ring, and 
an offset ring plus a pair of clumps. Some of these structures have slowly become revealed 
in numerical simulations of particular debris disks; we have chased them to their dens in the 
resonant landscape of the 3-body problem. The crude key we have assembled should help 
classify the debris disk structures seen by upcoming telescopes like SIRTF, SOFIA, ALMA, 
JWST and Darwin/TPF. 

Observing one of these structures instantaneously should allow us to categorize the 
planet as high or low mass (compared to Jupiter orbiting the Sun), and low or moderate 
eccentricity (compared to e ~ 0.2). In the case of a ring with a gap or an offset ring plus 
a pair of clumps, the image of the face-on cloud directly indicates the current location of 
the planet and points to its longitude of perihelion. In the case of a blobby eccentric ring, 
numerical modeling can potentially reveal these parameters. 

We thank Tommy Grav, Sean Moran and Mike Lecar for helpful discussions. This work 
was performed in part under contract with the Jet Propulsion Laboratory (JPL) through the 
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Table 1: Resonant Arguments. 
Resonance Nominal 1 /a e/ orce d/eo Argument Leading Term in F(a, e, eo) 
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"Includes contribution from indirect term. 



